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Abstract. For short-range interacting systems, no Schrodinger cat state can be stable when their environ- 
ment is in thermal equilibrium. We show, by studying a chain of two-level systems with nearest-neighbour 
Ising interactions, that this is possible when the surroundings consists of two heat reservoirs at different 
temperatures, or of a heat reservoir and a monochromatic field. The asymptotic state of the considered sys- 
tem can be a pure superposition of mesoscopically distinct states, the all-spin-up and all-spin-down states, 
at low temperatures. The main feature of our model leading to this result is the fact that the Hamiltonian 
of the chain and the dominant part of its coupling to the environment obey the same symmetry. 

PACS. 03.65.Yz Decoherence; open systems; quantum statistical methods - 03.65.Ud Entanglement and 
quantum nonlocality - 05.70.Ln Nonequilibrium and irreversible thermodynamics 
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1 Introduction 

The apparent classical behavior of macroscopic objects is 
thought to find its origin in the unavoidable interaction of 
any system with its environment [1,2]. For example, non- 
classical correlations between quantum systems [3,4] are 
expected to be fragile against this influence. This fragility 
of quantum entanglement is confirmed by studies show- 
ing that an initial entanglement between two independent 
open systems disappears in a finite time [5-7]. As a mat- 
ter of fact, when several independent systems interact with 
an environment in thermal equilibrium, their state generi- 
cally relaxes to their uncorrelated canonical thermal state. 
Consequently, even classical correlations are destroyed in 
this situation. This is not necessarily the case with a non- 
equilibrium surroundings. For an environment consisting 
simply of two heat reservoirs at different temperatures, 
the steady state of two non-interacting systems can be 
entangled [8]. Thus, in this case, an initially uncorrelated 
state can evolve into a quantum-mechanically entangled 
one. General results have been obtained concerning the 
steady state of a Markovian master equation of Lindblad 
form [9,10]. Within this approach, any pure state can be 
asymptotically reached with a purely dissipative dynam- 
ics, but not all states are attainable if the environment 
influence is required to be local. 

Probably the most amazing predictions of quantum 
theory arise when the superposition principle is applied 
to macroscopic objects. All states in the Hilbert space of 
any system can genuinely exist. There is no a priori re- 
striction, even in the case, for instance, of the Earth's cen- 
ter of mass. Quantum coherent superpositions of macro- 
scopically distinguishable states were first discussed by E. 
Schrodinger who considered a cat in a dead-alive state [11] . 



Here also, environmental degrees of freedom are thought 
to play an essential role. Under their influence, such freak 
states would decay very quickly into statistical mixtures if 
they would happen to occur [12-15]. Recently, Schrodinger 
cat states have been realized experimentally as superpo- 
sitions of motional wavepackets of trapped ions [16], mi- 
crowave cavity coherent fields [17,18], magnetic flux states 
of superconducting quantum interference devices [19], in- 
ternal states of trapped ions [20,21], photons polarizations 
[22], nuclear spins states in benzene molecules [23], free- 
propagating light coherent states [24], polarization and 
spatial modes of photons [25]. In all these experiments, 
the dissipative influence of the environment tends to de- 
stroy the created superposition of states and is one of the 
main obstacles to overcome to produce and observe it. 



However, considering the positive impact on quantum 
entanglement of driving the environment out of equilib- 
rium [8] , one can wonder whether a Schrodinger cat state 
can be stable in a multiple-heat-reservoir surroundings or 
in the presence of a monochromatic field. We address this 
issue in this paper by studying a chain of two-level systems 
(TLS), with nearest- neighbour Ising interactions, coupled 
to a heat bath and a monochromatic field, as illustrated in 
Fig. 1 . Couplings between the TLS are necessary to obtain 
a steady Schrodinger cat state [10]. It is proved below 
that, for any TLS system, its ground state cannot be a 
Schrodinger cat state when only short range interactions 
are present. Consequently, such a superposed state is nec- 
essarily unstable with an environment in thermal equilib- 
rium. We will see that, on the contrary, if the surroundings 
includes a monochromatic field, or a second heat bath, the 
TLS steady state can be a pure superposition of mesoscop- 
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ically distinct states, even if the interactions between the 
TLS are short-ranged. 

The paper is organized as follows. The model we con- 
sider is presented in the next section. Section 3 is devoted 
to the study of the TLS Hamiltonian. In the following 
section, we derive equations that determine the asymp- 
totic state of the TLS in the regime of weak coupling 
to both the heat reservoir and the monochromatic field. 
In section 5, we show that, when the field frequency is 
equal to a particular transition frequency of the TLS sys- 
tem, the TLS asymptotic state is, at low bath tempera- 
ture, a steady multipartite entangled state that cannot 
be obtained at thermal equilibrium. Moreover, a TLS- 
field coupling strength regime can exist, where it is a 
pure Schrodinger cat state. In the last section, we sum- 
marize our results and mention some questions raised by 
our work. The case of a two-heat-reservoir environment is 
considered in the last Appendix. 



2 Model 

We consider a system consisting of 2N two-level systems, 
a heat bath and a monochromatic field, described by the 
Hamiltonian 

H = H T ls + a^ 1 e f a^{a) + a) + ua)a 

2N 2N 

+ e CT > + £ «< + « + «) + H B , (1) 

n— 1 n— 1 

where u is the field frequency, e, e <C e and e/ « e are 
energies characterizing the coupling strengths of the TLS 
to their environment, a is a dimensionless constant that 
will be defined below, Hg is the bath Hamiltonian, tt and 
7r^ where v £ {x,y, z\, are observables of the bath, and 

27V-1 JV 

H T LS = - J ^2 a n a n+l ~ ^2 hn ( a n - a 2N+l-n) ( 2 ) 
n— 1 n— 1 

is the TLS Hamiltonian. The coupling constant J > and 
the fields h n are assumed to obey h n > and J2 n h-n < 
J/2. The annihilation operator a satisfies the bosonic com- 
mutation relation [a, a^] = 1. In the following, the nota- 
tion n is used for the eigenvalues of the number a^a. The 
Pauli operator has eigenvalues ±1 and the correspond- 
ing eigenstates are denoted by |±)„. The observables a* 
and are then defined by cr^|±)„ = |=p}„ and <jV 

■ ■ The main assumption underlying our model is that 
the dominant contribution to the interaction between the 
TLS and their environment is a uniform coupling to the 
heat reservoir, of the form ttS x where S x = J2 n a*- The 
term proportional to the energy e in (1), describes the 
small deviation, which is assumed local and generic, of 
the TLS-bath interaction from this ideal form. The cou- 
pling to the monochromatic field can be generalized to 
oT x tf J2 n v IT ,1(^ni' fl ' + Kiv a )- O ur m ain conclusions re- 
main the same provided Xnx~ i\N y — A/v+ix - i^N+iy 7^ 0. 




T 



Fig. 1. Schematic representation of a chain of two- level sys- 
tems with nearest-neighbour interactions and coupled to a heat 
reservoir of temperature T and to a monochromatic field of fre- 
quency LJ. 

The initial state of the total system is 

n = Y,r k i\k)(l\ ® Z- X e-**I T ® \a)(a\, (3) 
k,l 

where k and I run over the eigenstates of Htls, the den- 
sity matrix elements r^i are arbitrary, T is the tempera- 
ture of the heat bath, Z = Trexp(— H&/T), and \a) is a 
coherent state of the field, i.e., a\a) = a\a). Through- 
out this paper, units are used in which % = fee = 1< 
We assume that the average boson number a 2 is large. A 
time-dependent unitary transformation allows to change 
the Hamiltonian (1) into H + 2tf cos(ujt)afj where t is 
the time, and the initial state \a) into the field vacuum 
state. The time evolution of the TLS with these trans- 
formed Hamiltonian and initial field state, is identical to 
that ensuing from (1) and (3) [26]. The behavior we will 
find in the following, can thus be also obtained with a 
classical force in place of the quantum field. Another pos- 
sible physical interpretation of the large a limit, with 
the factor oT 1 in the Hamiltonian (I), is as follows. The 
strength of a local coupling between the TLS N and a 
cavity mode is proportional to V^ 1 ^ 2 where V is the vol- 
ume of the cavity, and hence, for a given energy den- 
sity waV" 1 , to a -1 . Thus, large a means large cavity 
with finite energy density [26]. The thermal average val- 
ues of the bath operators appearing in (I) are assumed 
to vanish, i.e., Tr(^7r) = Tr(J?7r£) = 0. This is the case, 
for example, for baths of spins or fcrmions in their high- 
temperature phases [27-30], or for heat reservoirs consist- 
ing of harmonic oscillators linearly coupled to the TLS 
[31]. In this last example, the bath Hamiltonian reads 
H& = J2q UJ q a \ a q where [a g ,ajj,] = 6 gg ', and the bath ob- 
servables 7r^ and 7r are linear combinations of the opera- 
tors a q and a q , and hence, the above mentioned average 
values clearly vanish. 



3 TLS Hamiltonian 

With the assumptions J > 0, h„ > and ^ n h n < J/2, 
_Eo = — J(2N — 1) is the ground energy of the TLS system 
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and its degeneracy is 2. The states 



it 



2N 2N 
® \+)n, |4)=® 



(4) 



constitute a basis of the corresponding eigenspace. More 
generally, the configuration states \rj) — ® n \r] n )n where 
Vn G {+!~ } and T] = (771 , . . . , t)2n), are eigenstates of 
Htls- An important property of this Hamiltonian is that 
it commutes with the symmetry operator II defined by 



II|?7) = ®„ I r? 2 jv+i-n )n> 



(5) 



where ~EE~ = =p. For example, for N = 2, II | + + H — ) = 
I H ) and XI | + H ) = | + H ), with a simpli- 
fied notation. Clearly, II 2 = 1 and II has only two eigen- 
values, 1 and — 1. The Schrodinger cat states 



|Scs-)=2- 1 /2(|^) T |4)) 



(6) 



are ground states of Htls and eigenstates of II. For N — 2 
and general values of the coupling constant J and fields 
h n , Htls has ten eigenvalues E k . Four are nondegen- 
erate and correspond to configuration states which are 
eigenstates of II. The six other ones are doubly degen- 
erate. For larger values of N, there are levels with higher 

degeneracies. For example, for N = 3, | — I ), 

I - + + H ), I + + + H h) and | -M h) are 

eigenstates of Htls with the same eigenenergy — J— 2h-2- 
Specific energy levels will play an important role in 
the following. Consider a configuration 7] with a single in- 
terface, i.e., such that rj n — r\ for n smaller than a given 
integer m < 2N, and rj n = rj for n > m. The correspond- 
ing eigenergy is Ek = Eo+2J—2i]J2™ = i h n where m! = m 
if m < N and m! — 2N — m otherwise. For general values 
of J and h n , the degeneracy of such a level is, for any N, 
1 if m = iV and 2 otherwise. These levels are the lowest 
ones, see Fig. 2. In this figure and in the following, we use 
the notations 



Zfcl 



IN+l 

-)n. 



■,\b) = \-h 

M ± )=i±\ 



1 N- 



-i.-.,(7) 
l-)iv---, 



where I± = 2 _1 / 2 (1 ± II) and . . . means that the state of 
TLS n is that of TLS n — 1. These states are eigenstates 
of both Htls and II. The corresponding eigenenergies are 
denoted by E a , E bl E c and E d . 



4 TLS asymptotic state 

For a large average boson number a 2 , there exists a time 
regime where its variation is negligible and the TLS chain 
reaches an asymptotic state poo [26] . To determine poa , we 
first write the reduced density matrix of the system made 
up of the TLS and the monochromatic field, at positive 
times t, as 



Ptls+fH) - 2~ 



dze- lzt Tir- 



(z-cy 1 n 



(8) 



Ed 



+ + >. 



-+++>+! — +> 



-+++>- 



+> 



7 1 



+ — >+|+++-> J+ — >-!+++-> 



++++>+ 



> 



++++>- 



> 



Fig. 2. Lowest lying eigenstates of both the TLS Hamiltonian 
Htls and the symmetry operator n for 4 two-level systems. 
The eigenenergies are Eo = —3 J, E a ~ — J — 1h\ — 2fi2, E c — 
-J-2h,!, E d = -J + 2/ii and E b = - J + 2hi +2h 2 . The states 
on the left (right) correspond to the eigenvalue 1 (— 1) of IL 
The arrowed lines connect the 'left' and 'right' states which 
are coupled by the monochromatic field. 



where Trg denotes the partial trace over the heat bath, 
£ is a positive real number, and the Liouvillian C is de- 
fined by C ... = [H, . . .]. The matrix elements fkhis(z) = 
(k\(h\Tris[(z — C)^ 1 (2]\l)\S) of the Laplace transform of 
Ptls+f, can be written as 



^ knl 



Is.k'n'l' s 



S'(z)r k >, 



(9) 



where the functions rknis,k'n'i's'(z) depend only on the 
heat bath part of the initial state (3) [8]. The right hand 
side of this equation can be read as the product of a square 
matrix T(z) with a column vector. A master equation for 
Ptls+f can then be derived with the help of the inverse 
matrix £ = see Appendix A. 

We are concerned with the limit of weak coupling of 
the TLS to their environment and with e, e/ <C e. It is only 
in the limit of weak coupling to its surroundings, that the 
density matrix of an open system with Hamiltonian H., re- 
laxes to the canonical equilibrium state p eq cx exp(— %/T) 
when its environment is in equilibrium with temperature 
T . To obtain in the regime of interest, the expansion 
of S to first order in e/ and to second order in e and e, is 



required. We find S 
where 



+ a- 1 e f V TLSF + . 



khls.k'fi 

yTLS B 
khls 7 k' ft 

yTLS F 
khls.k'fi 



= [Z - U} ki + (S - h)u>}5 k i k 5lil5n>n5 Sl s (10) 

= -ijki,kn>(z + (s-n)ui)5fi'n.Ss's (H) 

= Sk'k<7l'l8h'ny/s'8s's+l + sSg's-1 (12) 



— Sl'lVkk'fis's^/nfiri'fi-l + n'&n'h+l- 

In these expressions, the notations u k i = E k — Ei and 
Cfci = {k\a%\l) have been used. The functions Jkl,k'l', 
which describe the influence of the heat bath, can be ex- 
pressed in terms of bath time-dependent correlation func- 
tions, see Appendix B [8,32,26]. For a large number a 2 
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and a finite time t, the matrix element (k\(h\prLS+F\l)\s) 
is non- negligible only if n and s are close to a 2 . Conse- 
quently, these integers can be replaced by a 2 in (12). Since 
they appear in (10) and (11) only via the difference s — n, 
it is convenient to define 

f u( Z )=^Z f knln+p{z). (13) 

n 

For p = 0, this expression gives the matrix elements of the 
Laplace transform of the TLS chain state 

p(t) = Tr FPTLS+F (t) (14) 

where Tip denotes the partial trace over the monochro- 
matic field. 

As we are interested in the TLS asymptotic state, 
we write the column vector r(z) of elements fknis{z) as 
r(z) = ^2 x r x /(z — x) + r(z) where x is real and r(z) 
has no singularity on the real axis. The potential pole at 
z = x, gives an undamped component of the TLS state, 
of frequency x, and hence contributes to poo- Since the 
elements of the vector £(z)r(z) are equal to matrix el- 
ements of the TLS+field initial state, it is non-singular, 
and thus the vector must fulfill S(i0 + + x)r x = 0. Let 
us focus on the term x = 0, which always exists, and de- 
note by u the zeroth order of r in an expansion in terms 

of e, e and e/. Similarly to (13), we define uff from the 

elements of u. From S°(i0 + )u = 0, it ensues that uff is 
non- vanishing only if u>ki = poJ- Then, from the next order 
terms of S(z0 + )r = 0, we obtain 

Er / (p+i) , (p-i)\ t (p+i) , cp-i)\' 
Kj + u kj ) - Vkj (u), + u}, ') 

3 

- i IW (* 0+ + P") u i% = 0- (15) 

k',V 

where k, I and p are such that wm — pw and u k ,/, = 
if oJk'V 7^ p'u. If uj is chosen arbitrarily and is not equal 
to any ojki jp then the only non- vanishing u k ^ obey p = 
and Ek = Ei, and arc solutions of equation (15) with 
no monochromatic field term. In this case, the field influ- 
ence manifests itself only through higher order corrections. 
In the opposite case, for a resonant field frequency, the 
dominant contribution to the asymptotic behavior is af- 
fected by the coupling to the monochromatic field. For the 
more general TLS-ficld coupling a _1 e/ J2 n v a ni^nvo) + 
Xn V a), in the first line of equation (15), Uki is replaced 
by J2 n v ^nv(k\<Tn\l) in the first and third terms, and by 
J2n v Kiv{H a n\l) m the second and fourth ones. 

Since iiff = if uj^i ^ 0, the matrix elements (k\p\l) 
of the TLS state (14) such that Ek ^ Ei, have no steady 
component. Moreover, it is clear from the expression (10), 
that, for uoki = 0, the only possible undamped component 
of {k\p\l} is of zero frequency, and hence 

{k\ Poo \l)=u k y (16) 



for Ek = Ei. Actually, the other non- vanishing u k p / give 
also contributions to p^. This can be seen as follows. Let 
us define the vector rp w of elements (r ) kms+p- It satisfies 
S(i0 + + Pbj)rpu = with the approximate £ discussed 
above, and thus, corresponds to a pole of r(z) on the real 
axis at z = Puj. Moreover, Y,h( r P^)kAih = Es( r o)feSift+P 
is equal to u kl to zeroth order, and hence (16) gener- 
alizes to (fc|/Joo|0 = u k P ^ exp(-ipujt) for Ek = Ei + puj. 
In Appendix A, we derive from (10)-(12), a Markovian 
master equation for p(t) which reduces to Redfield equa- 
tion when Ef = 0, and show that a time-periodic p = 

J2 P ex P( — ipojt)g^ satisfies this equation only if = 

u k P i to lowest order in e, e and e/. In this Appendix, we 
also show that, in the case of a single TLS coupled to 
a zero-temperature bath, (15) gives the same asymptotic 
state as that obtained from the optical Bloch equations in 
the rotating wave approximation [26] . 



5 Schrodinger cat regime 

In this section, we consider the case of a field frequency 

uj = 2h N = E c -E a = E b - E d , (17) 

sec Fig. 2, and of a low heat bath temperature, 

T < h n < J. (18) 

We will see that, in this limit, the asymptotic state p,^ is 
a steady state of the form 

Poo =p|Scs^)(Scs-| + (l-p)|Scs + )(Scs+| (19) 

where the Schrodinger cat states |Scs T ) are given by (6). 
For vanishing coupling to the monochromatic field, e/ = 0, 
p = 1/2 and Poo = (| ft) (ft | + | ft)(ft|)/2 is simply the 
zero temperature thermal state corresponding to the TLS 
Hamiltonian Htls- Fot Cf 3> e, p — 1 and the TLS asymp- 
totic state is a pure Schrodinger cat state. At the end 
of this section, it is shown that, provided p ^ 1/2, the 
mixed state (19) is not as highly entangled but remains 
nevertheless multipartite entangled. To obtain (19), we 

first define p 8 = ^2 k l u k °i\k)(l\, which is the steady com- 
ponent of poo. From (16) and Cauchy-Schwarz inequal- 
ity, it follows that any matrix element of poo satisfies 
KfclPoolOI 2 < (k\ps\k)(l\p s \l), and that (k\ Poo \l) = (k\p s \l) 
if Ek = Ei. In the following, we focus on p s and show 
that its only nonvanishing elements in the limit (18), are 
(Scs =F |p s |Scs =F ), which leads to the result (19). 



5.1 No thermal equilibrium Schrodinger cat 

Before embarking on the determination of poo, let us first 
show that the thermal equilibrium state p eq of any TLS 
system cannot be of the form (19) with p ^ 1/2, in the ab- 
sence of long range interactions. At temperature T, p eq cx 
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J2k \k)(k\ exp(— Ek/T) where |fc) and Ek are the eigen- 
states and eigenenergies of the TLS Hamiltonian H. If the 
temperature is too high, p eq is a mixture of a large num- 
ber of states. At zero temperature, p eq is the equal-weight 
mixture of the ground states of H. Thus, in this case, p eq 
can be of the form (19) but with p = 1/2 only. The ther- 
mal state is given by p eq = p|0)(0| + (1 — p)|l)(l| with 
p 7^ 1/2, only if |0) and |1) are the ground state and first 
excited state of H, and are non-degenerate, and only in 
the limit E 2 3> T. We thus assume that H has a nonde- 
generate ground state |0). When no long range interaction 
is present, H can be decomposed as H = Hi+H' where no 
observable of TLS 1 appears in H' and there exists a TLS 
n 7^ 1 which is not affected by Hi- It can be proved that 
|0) cannot be |Scs T ) by evaluating the average value A = 
(Scs^-HIScs*) as follows. Since (WIH) « i< - |+>i = 0, 
and similarly for Hi, A = [(^\H\ty) + (fr \H\ fr)]/2. Conse- 
quently, A > (0|"H|0), as this last value is the minimum 
possible one for the average energy, and |0) cannot be both 
Iff) and I-U-). Thus, |0) is not |Scs T ). This result generalizes 
to all Schmidt decomposable states, see Appendix C. 

5.2 Uniform coupling to the heat bath 

As we are interested in the limits e, e/ <C e, we first con- 
sider the case e = Cf = 0. For these particular values, 
the total Hamiltonian (1) commutes with the symmetry 
operator II defined by (5), and there thus exist several 
TLS asymptotic states. This can be seen as follows. Con- 
sider that the TLS initial state ^ fc l rki\k)(l\ is a statis- 
tical mixture of eigenstates of LT with eigenvalue 1. Since 
[H, LT] = 0, the TLS reduced state at any time is also 
such a mixture. The conclusion is obviously similar for 
the eigenvalue —1 and hence the asymptotic state cannot 
be the same for these two kinds of initial states. Note that 
this argument is valid for any coupling strength e. We 
remark that, though Htls and S x = Yln a n both com- 
mute with the symmetry operator II, they do not have, 
for general values of J and h n , any common eigenvector, 
and hence that there is no decoherence-free subspace in 
the TLS Hilbert space [33], see Appendix D. 

For ef — 0, equations (15) separate into the inde- 
pendent equation sets J2 k , v jki,k'i'(i0 + + P L0 ) u k'l' = 0) 

(p) 

where k and I are such that lum = pw, and u k ,\, = if 
uJk'V 7^ pu- The equation set p = gives the steady com- 
ponent p s = l u. a )\k)(l\ of poo. Up to now, the basis 
set {|fc)} has not been fully specified. The states \k) are 
eigenstates of Htls, but, for a degenerate energy level 
E k , infinitely many choices are possible. To obtain the ex- 
pression of p s , we consider the diagonalized form of this 
density matrix in the appropriate unknown basis set {\k)}. 
The equations 

£7«,k'k'(iO + )uj&,=0 (20) 
k' 

where k and I are such that Ek = Ei, then determine both 
{\k}} and the populations u k ° k . Let us introduce the coeffi- 
cients T k i = jkk,u( i0+ ) cxp(-Ei/T). They obey T ik = T kl 



and Y^k ~ 0> see Appendix B, and hence, for any set 

Wk}, 

^2r k npk(pi = -2^2 T ki(<fk - <pi) 2 - (21) 

k,l k,l 

By definition of Y k i, this sum vanishes for (p k = u k ° k x 
exp(E k /T). Since T ik = T u and T u > for k ^ I, the 
vanishing of (21) is equivalent to TkHPi = 1 . As shown 
in Appendix B, for k ^ I, Tki is proportionnal to Ski = 
(k\S x \l) 2 when e = 0. Since the right hand side of (21) is 

equal to zero only if ipk = <Pz when Tki ^ 0, u k ° k /u^ = 
exp(— LOki/T) when Ski ^ 0. This relation ensures that the 
equations (20) with k ^ I, are satisfied, see Appendix B. 

Thus, p s is determined by the equations (20) with k — 
I, or, equivalently, by the vanishing of the sum (21). It can 
always be written as 

q q ke£ q 

where z q = J2kes q exp(-E k /T) and J2 q Pq = 1- This ex- 
pression satisfies (20) if the states \k) and the sets £ q are 
such that Ski = for any k e £ q and I e £ q i where 
q' q. Obviously, as soon as there exist more than one 
set £ q , the decomposition (22) is not unique, as a new 
family of sets £ q can be simply defined by unioning sets. 
Thus, we assume that a set £ q cannot be divided into sub- 
sets £ q i and £ q » such that Ski = for any k € £' q and 
I € £ q ". In this case, every state \k) in (22), is eigenstate 
of n and the value of the corresponding eigenvalue de- 
pends only on the set £ q to which k belongs. This can be 
shown as follows. Consider the subspace S q of the TLS 
Hilbert space, spanned by {\k}} where k £ £ q . It can al- 
ways be decomposed as S q = S q+ ® S q - where S q ± is 
spanned by the states (1 ± n)|fc), which are eigenstates 
of both Htls and n. Since [S X ,U] = 0, (^\S X \^') = 
for any \i[>) e S q+ and e S q -. As £ q cannot be sub- 
divided, as assumed above, S q = S q+ or S q = S q -. In 
other words, all the states \k) given by £ q , are such that 
II|fc) = ±|fc). For N = 2, the number of levels Ek is small 
enough to allow the complete determination of the min- 
imal sets £ q . For general values of J and h n , there are 
only two sets, corresponding to the two eigenvalues of LT. 
For hi = 0, there are three sets, one giving only the state 
(| + - + +)-|-- + -)-|+ + -+) + |- + --))/2.We 
remark that, since this state is eigenstate of both Htls 
and S x , it is decoherence-free, see Appendix D. The situ- 
ation is similar for hi — and hi = hi. 

For our purpose, we only need to know a few states \k) 
and the corresponding two sets £ q . The two ground states 
|Scs~) and |Scs + ) are eigenstates of n with eigenvalues — 1 
and 1, respectively. Thus, they are given by two different 
sets, £- and £ + . Consider the doubly degenerate energy 
levels E e and Ef corresponding, respectively, to the states 

1 The matrix of elements T k i is real and symmetric and hence 
diagonalizable. Using (21) and Tu > for k I, it can be 
shown that its eigenvalues are negative, which leads to the 
equivalence. 
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i|+) 2 • • ., where 
± 



|e±> =X±|+> 1 |-) 2 ... and |/± 

X± = 2- 1 /2(i ± n) 2 . Since 5 Scs±e± = 5 Scs±/± = 1, e 
and / belong to £±. By evaluating the matrix elements 
of 5a; between the one-interface states discussed at the 
end of section 3, and noting that the degeneracy of a one- 
interface level is 1 or 2, it can be shown that £+ gives all 
the states l+\fj) where \rj) is a one- interface configuration 
state which is not eigenstate of IT, \a), \b), and |Scs + ). 
Similarly, £_ gives all the states I-\rf) and |Scs _ ). 



5.3 General coupling to the heat bath 



The non- uniqueness of p s for e = e/ = stems from 
the particular form of the coupling to the environment 
in this case. If the uniformity of this coupling is bro- 
ken, even slightly, this indeterminacy disappears, as we 
show here. For e / and e/ =0, the potential sets £ q , 
appearing in the decomposition (22), are determined by 
e(k\S x \l)(A\n\B)+ej: n JkK\l)(AK\B) = where \A) 
and \B) are any eigenstates of the bath Hamiltonian 
since T^i vanishes only if this condition is satisfied, see 
Appendix B. In general, it is equivalent to (k\a!^\l) = 
for any n and v. In other words, k and I belong to the 
same set £ q if (k\a^\l) ^ for some n and v. As a con- 
sequence, k and / belong to the same set £ q if there exist 
Ji, J s , ni, n s , n, v\, v s and v such that 
® r (c^|>)(Jr|)CnlO 7^ 0' which means that there is a 
bath- induced transition path between the states |fc) and 
\l). Since 



r=2 



E I^KIiOllOV-iKI^O.KIO 

(k\K=i«\i) , (23) 



> 



k and I belong to the same set £ g if there exist n\, . . . , n s , 
n, v\, . . . , u s and v such that (fc|<g> r (T^(7^|/) ^ 0. Consider 
any state \k) appearing in (22). It can be expanded on the 
basis of configuration states \ff ) — ®„ \rj rn ) n as |fc) = 
^2 r K\l] r )i which can be rewritten as \k) = J2r ®n 
( a n) Tr "\ a ) where \a) = \+)i ■ ■ ■ \— )jv+i • • • , and r rn is equal 
to 1/2 + r\ rn j2 for n > N and to 1/2 — rj rn /2 other- 
wise. As the eigenenergy E a is nondegenerate, \a) is nec- 
essarily present in (22). Since at least one A r is nonzero 
and (a\ <8> n (a^) Trn \k) = X r for any r, a and k belong 
to the same set £ q . Thus, there is a unique set £ q and 
ps = peq oc exp(—Ek/T)\k)(k\ is the thermal state 
of the TLS system. Consequently, in the low temperature 
limit (18), Poo = Ps = (|fr)(fr| ± I TO DA which can be 
written as (19) with p = 1/2. We remark that the above 
proof applies even if = ix z n = 0, and that the only re- 
quired assumption on the parameters oIHtls is that they 
are generic enough that E a is non-degenerate. 



5.4 Low temperature limit 

In this section, we show that, for any ef and any e ^ 0, 
(k\p s \l) vanishes in the low temperature limit (18), if Ek = 
Ei is larger than E + 3 J. This result, together with those 
of section 5.2, lead, for e, e/ <C e, to the form (19) for p^. 

Let us first write p s = ^2 k u< kk \k)(k\ in diagonal form, and 
observe that, if E\ < Ek, 7fcfc,ii(«0 + ) vanishes at low tem- 
peratures, see expression (40), and that, for any I such 
that Ei > E a , there exists k such that Ek < Ei and 
7fcfc,H(i0 + ) ^ 0. In other words, the upward transition 
rates vanish at zero temperature, and there is a down- 
ward transition from any excited state but the first one 
| a) = |+)i . . . |— ) n+1 ■ ■ ■■ This last property can be seen as 
follows. Except the ground states | ff) and | JJ-), and the one- 
interface states |±)i|=f) 2 ... and ■ ■ .\±)2n-i\t)2N, any 
configuration state \rf) presents a sequence |±) n |=F) n+ i 
where n </ {l,2N - 1}. Thus, if \rj) is neither one of 
the states mentioned above, nor |±)i . . . |=f) n + i . . ., it con- 
tains a sequence | = F)n-i|i)n|T)n+i where n </ {l,2iV — 
1}, or a sequence |±)„„i|±)„|=f)„ + i|=f)„ +2 where n # 
{1,N,2N — 1}. In the first case, the energy of a^lV) i s 
lower than that of \rj) of —2J±2h n where h n = —fi2N+i-n 
forn > N and h n otherwise. In the second case, the differ- 
ence between the energies of \rf) and of er^|7?) is ±2h n , and 
that between the energies of \rj) and of cr* +1 \r]) is ^f2h n +\. 
Consequently, since h n > 0, there always exists a configu- 
ration state (TnlV) with one TLS flipped, of energy lower 
than that of \vj). This result can be extended to the states 
|±)i|=F>2 • • • and . . . |±) 2 jv-i|t)2JV, since flipping, respec- 
tively, the first and the last TLS, leads to one of the two 
ground states |f|-) or |JJ.), and to \b) = . . . |+)jv+i • • 
since the energy of a^\b) and cr^ +1 |6) is Ed = Eb — 2Hn, 
see Fig. 2. Thus, \a) is the only exception. Consider now 
any / ^ a and expand the corresponding state on the ap- 
propriate configuration state basis as \l) = Yl^—i ^r\f] r )- 
As seen above, there exists n such that the energy Ek of 
<Jn\f]\) is lower than E\. There are possibly s' states cr^|^ r ) 

of energy E k . The state \ip) = X)r=i ^nlVr) belongs to 
the TLS Hilbert subspace spanned by the states |fc) of en- 
ergy Ek- Consequently, there is at least one |fc) such that 
(k\a^\l) — (k\ip) 7^ 0, and hence, such that 7fc/ Cj o('i0 + ) 7^ 0. 

We now consider equation (15) with p = and T = 0. 
Any eigenenergy E k is of the form Ek = Eq + 2mJ + 
Y] h n T n where r n G {— 2, 0, 2} and m is the number of 
interfaces of the corresponding configuration states. Con- 
sequently, for general values of J and h n , if ujki = ±w = 
±2/i jv then k and I correspond to the same number m 
of interfaces. In particular, the monochromatic field does 
not couple one-interface states to multi-interface states. 
Therefore, at T = 0, 



{ a 'k(. U kl 



(i) 



ieA r 



* E ^ u u ] - l^tu) + * E WkH - 0, (24) 



ieA L 



For TV = 2, e = c and / = d, see (7). 



where jki is the zero temperature limit of 7fefe,i;(«0 + ), Ah = 
{k : E k > E Q + 3J}, A L = {k : E k < E Q + 3J}, and 
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k € Ah- The solution to this equation set can be in- 
terpreted as the asymptotic state of a system evolving 
under the influence of a zero temperature heat bath, a 
monochromatic field, and an additional decay mechanism 
characterized by the rates ^2 ley ^ L Jik- Equation (24) leads 

t° J2keA H .ieA L Hk u k°k = 0- Since all the coefficients ji k in 
this sum are positive, u k ° k = if one ji k is different from 
zero. The results obtained above show that this is the case 
for k G Ah corresponding to the lowest E k given by this 



,(o) _ (±1) 
3 ' a kk ~ u kl 



u\t^ = for these k, and, in the 

(0) 



'■Ik 

equations (24) determining the other u k " k ' where k e Ah 
the sets Ah and Al can be replaced by the sets A' H and 
A' L obtained by, respectively, removing these k from Ah 
and adding them to Al- This new equation set leads to 

u k°k = u< id^ = u lk^ = for A: e A' H corresponding to the 
lowest Ek given by this set. By repeating this procedure, 

it can be shown that u k ° k = for all k e Ah, i.e., such 
that E k > E + 3 J. 

For e, €f <C e, we know that p s is given by (22). The 
results of this section show that the probabilities p q where 
q 7^ ±, vanish at low temperatures, since the correspond- 
ing sets E q give states with more than one interface. Thus, 
the low temperature TLS asymptotic state is of the form 
(19). By slightly modifying the above derivation, it can 
be shown that p q <C exp(— 3J/T) where q ^ ±. At low 
temperatures, jkk,u( iQ+ ) — Hk exp(-u kl /T) if Ei < E k , 
see expression (40). Hence, the terms neglected in equa- 
tion (24), are approximatively given by — *X^e.4 L 7'fc x 

exp(-w fe/ /T)4 0) . Since u<°> ~ p ± exp[-(Ei - E )/T] for 
! £ f±, and Al C £+ U the above sum is far smaller 
than exp(— 3J/T), and hence exp(3J/T)u k ° k where k e 
Ah, vanishes in the zero temperature limit. This result 
will be useful in the following. 



5.5 Steady Schrodinger cat state 

We assume here that t < tj < e. To obtain p s in this 
case, it is convenient to first rewrite equations (15) in ma- 
trix form as G W uW + e / G^ ) u^ +1 ) + e f G^lu^-^ = 
where is the column vector whose elements are the 
non-vanishing u kl , i.e., such that u k i = pu. In the follow- 
ing, we use the inverse matrices (Gq^) -1 where p / 0. 

Their existence can be shown as follows. If were an 

(v) 

eigenvalue of G , possible asymptotic states for e/ = 0, 
would be poo = p eq + xp osc where 

Pose is a matrix with 
only nondiagonal elements, oscillating at frequency pu, 

constructed from an eigenvector of G ' with eigenvalue 
0, and x is any real number. But, since p^ is a density 
matrix, its off-diagonal elements obey x 2 |(/c|p osc |l)| 2 < 1, 
which is not possible for any x if p osc ^ 0. Thus, is not an 
eigenvalue of Gq*\ which is hence invertible. Physically, 
this means that the possible oscillating components of p^ 
with frequencies pu, are induced by the monochromatic 
field and disappear when the coupling to it vanishes. As 
a consequence, the contribution to p^ which oscillates at 



frequency u has an amplitude of the order of e/, and, more 
generally, u( p ) is of the order of e^ 3 '. Thus, u^°) satisfies 



e)G 



i 0) {< 



G 



(i) 



-e)G 



(0) 



(-i) 



u<°> = 0, (25) 



where terms of order and higher have been neglected. 

We denote by Gg , the matrix G with e set to zero. 
As seen in section 5.2, Gg ' has several eigenvectors 
with eigenvalue 0, whose elements are z~ l exp(—E k /T)5hi 
if k e £ q and otherwise, with the basis set {\k}} con- 
sidered in section 5.2. There also exist column vectors <I> 9 

which obey ^^Gg ' = and <fr^'J' 9 < = 5 qq >, and whose 
elements are Ski if k S £ q and otherwise, see Appendix 
B. To zcroth order in e/ and e, u^ ' = J2 q Pq^q- Replac- 
ing u(°) by this expression in (25) and neglecting terms 
proportional to powers of e, leads to Y] , r qq rp q r = with 



' qq> — ^ q 



G (o) F (i) G W + G (o) F (-i) G ( 



(-i) 



(26) 



where F^ 1 * 1 = (G u ±1 '') _1 . Since the monochromatic field 
couples states given by different sets £ q , sec Fig. 2, the 
probabilities p q appearing in (22) are no longer indepen- 
dent from each other. 

With the explicit expressions of the matrices G±\ and 

G^ 1 \ and of the vectors <& q and ^ q , we find 



' qq 



-E„,/T 



79' 



ds:k'] (p) (p) F (p) 
-°q d q'\ a kl a k'l' b , 



k.l,k' ,!' 
= ±1 



(27) 

where o^f^ = (k\a^\l) when u k i = ±u and otherwise, 
and 5 q = 1 when k € £ q and otherwise. For the more 
general TLS-field coupling oT x e s Ylmv a n(^n^ + A*„a), 



t (1) 

7 kl 
J 1 ) 



.(-1) 



( a ik )*! which gives 



T, n ,v X nv(k\<\l) and a kl 
= 2~ 1/2 (\ Nx - i\ Ny - X N+ i x - i\N+i y )- To obtain 



the expression (27), we have used the fact that F^ fe/; , 
vanishes if k and k! do not belong to the same set £ q , 
or if I and V are elements of different sets, and Fjj.7fe'J' = 

0?lki'k')*> see Appendix B. The coefficients (27) satisfy 
J2 q r qq' = 0. We have seen in the previous section that p q 
where q ^ ±, vanishes at low temperatures, faster than 
exp(— 3J/T). The populations p- and p + obey r__p_ + 
r_ + p + + ^2 q jL± r -qPq = 0. For low temperatures, r± + ~ 
± exp(— E a /T)r and r T _ ~ ± exp(— E c /T)r where 



r = e Eo/T Re 



(lc-a,c-a(-i0 + +W)) 



(28) 



since F, 1 ? _ cx S kc -Si a at zero temperature, see Appendix 
B. Consequently, ~ — r jr |_ ~ cxp(— u/T) van- 

ishes in the low temperature limit (18), and thus the asymp- 
totic state poo = |Scs~)(Scs - | is a pure Schrodinger cat 
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state. We remark that, at strictly zero temperature, there 
is no upward transition depopulating the ground state 
|Scs + ) and hence the TLS system does not necessarily re- 
lax into the pure state |Scs~). A steady pure Schrodinger 
cat state can also be obtained by coupling the TLS to a 
second heat bath instead of a monochromatic field, see 
Appendix E. 

5.6 Steady multipartite entangled state 

As shown in the previous section, the TLS asymptotic 
state poo is a pure Schrodinger cat state when e/e/ <C 1. 
For finite values of this ratio, p^ is given by (19), and 
is not such a pure superposition of mesoscopically dis- 
tinct states. However, it remains multipartite entangled 
as long as p ^ 1/2. More precisely, there is no partition 
of the TLS system, with respect to which, it is separable. 
It is enough to prove it for an arbitrary bipartite split- 
ting. Let us then consider such a partition, i.e., any two 
subsets of {1, ... , 2A}, and name TLS\ and TLS2 the cor- 
responding TLS systems. The asymptotic state of TLS\ is 
pi = TrTLS 2 /°oo, and that of TLS2 is obtained by tracing 
over TLS\. The expression (19) can be rewritten as 

poo = + iwi) + {\-p) (iwi + iwi). (29) 

which shows that p T is an equal- weight mixture of the all- 
spin-up and all-spin-down states of TLS T . Consequently, 
the von Neumann entropy of p T is S(p T ) = — Tr(p T In p T ) = 
In 2. On the other hand, the entropy of the complete TLS 
system state (19) is S(poo) = —p\np—(l—p) ln(l— p). For 
p 7^ 1/2, S(p T ) > S(poo), and hence, the systems TLS\ 
and TLS2 are entangled [34] . 

6 Conclusion 

In this paper, we have seen that, for any TLS system 
with short-range interactions, no Schrodinger cat state 
can be stable when the system environment is in ther- 
mal equilibrium. To examine whether this is possible when 
the environment is out of equilibrium, we have studied a 
chain of two-level systems coupled to a heat reservoir and 
to a monochromatic field. For any even number of TLS, 
we found a regime of Hamiltonian parameters where the 
asymptotic state of the TLS chain is a pure Schrodinger 
cat state at low temperatures. Though obtained for a 
specific model, this result, together with that of Ref. [8], 
suggests that, more generally, driving the environment 
out of equilibrium can enhance considerably non-classical 
features of an open system. It would be interesting, for 
other models, to study how diverse non-classicality crite- 
ria, based on Wigner function for example [15], or entan- 
glement measures, evaluated for the system steady state, 
change with the distance from equilibrium of the environ- 
ment. 

The existence of a Schrodinger cat regime may not 
be specific to the TLS chain studied in this paper. It en- 
sues from some main features other systems can present, 



which are the following. The ground level of the system 
is degenerate. This is essential since, as we have seen, a 
non-degenerate ground state cannot be a superposition 
of macroscopically distinct states. The system Hamilto- 
nian and the dominant component of the interaction with 
the environment commute with the same symmetry op- 
erator. This leads to the possibility of unequal steady- 
state populations for equal-energy eigenstates of the sys- 
tem Hamiltonian. The relation between the populations 
of the ground states is determined by the non-symmetric 
part of the system-environment coupling. The presence of 
non-degenerate energy levels plays an important role in 
the existence of a regime where one of these two popu- 
lations vanishes. It would be of interest to study other 
systems showing the same features. 



A Markovian master equation 

By inverting the matrix T(z) in (9), one obtains the master 
equation 

^ Q/^ + S 2 

Ekhls,k'n'VS'{z)rk'n'i's'{z) = r k i—==e~~ a (30) 
,, frf. vnlsl 

k',n',l ,s' 

where S = r _1 . Expanding this matrix as £ = S° + 

-ZTLSB + a -l ef STLSF + where ^TLS B and 

S TLS F are given by the expressions (10)-(12) and tak- 
ing into account that the number a 2 is large, leads to 

d t r { k f + i( Wfcl - pu)r { $ 

- V f dt'Kuvv (t')e ipujt 'r { k % (t -t') = (31) 
k >,i> Jo 

where r[f = Ysn(M(n\pTLS+F\l)\n +p), which is related 
to (13) by Laplace transform. The time functions Kki.k'i' 
are given by 

K khk , v {t) = e u ^C l k l(t) + e^'CHi-t) 

- E [5u> e lt ^ C% (t) + 5 kk , (-*)] (32) 

3 

where C k lii(t) = Tr{f2ir k i (t)n k 'i>] are bath correlation 
functions [26,8], and, with the Hamiltonian (1), 

n kl =eJ2(k\<\l)7r + eJ2(kK\l)<, (33) 

n n,u 

and ir k i(t) = exp(— UH^-km cxp(itiJg). It is clear from 

the above definition that the components r$ = (k\p\l) 
are the matrix elements of the TLS reduced density matrix 
P = T^fPtls+f- 

For a large boson number a 2 , the time evolution of 
the monochromatic field is essentially not affected by its 
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interaction with the TLS and rff ~ exp(ipwi)rj^' 1 . More- 
over, since the coupling of the TLS to the heat reser- 
voir is weak, the Markovian approximation rffl (t — t') ~ 

exp(iuikit')r^ (t) can be used in the right hand side of 
(31). These two approximations give the Markovian mas- 
ter equation 



,(°) 



If the correlation functions C k ; t , vanish fast enough in the 
long time limit, the functions (37) are non-singular on the 
real axis where they can be written as 



d t r kl + ioj k ir k i - ^ -fu,k'i' (i0 + + u k >v Vk'V 

k',l' 

- 2ie f cos(wi) ^ {ojlTkj - ^kjrji) = (34) 

3 

where r k i — r k °i = (k\p\l) and the functions jhl,k'l' are 
given by (37), which reduces to Redfield equation [32,26] 

for e/ = 0. Assuming time periodic r k i = J2 P cxp(-ipujt) g k ^ 
leads to 



lki,k'i'(i0 + + cj) = 7T 

A,B ^ 

x S(E A - E B + oj + u> Vk ) + 5(-E A + E B +uj + uji k ,) 

- 5 W Y,( A MB){B\Tr jk ,\A)5(E A -E b +lj + w w ) 
j 

-5 kk ^{A\^ Vj \B){B\TT J i\A)5{-E A + E B +oj + Lo jk )^, 



i(cv k i - pu)g k p { - 7w,fe'/'(*0 + + ^k'l')Q k % 

k',V 

E\ ( (P+l) , (P-1)\ / (P+1) , (P-1)\ 

[<Wfej + Qkj ')-^kj{Q)i '+Q)i ') 

3 



(35) 



for any k, I and p. We solve this equation set pertur- 
batively in both the coupling to the heat bath and to 
the monochromatic field, and name uff the zeroth order 



of ■ The first term of (35) imposes that ujy 
uj k i 7^ pui. Then the lowest order of (35) for k, I and p 
such that uj k i = puj, gives equation (15). 

For a single TLS coupled to a zero-temperature heat 
bath and to a field of frequency uo = 2h, described by the 
Hamiltonian H = — ha z + a^ 1 ef(7 x (a J< + a) +uja J< a + 7ra x + 
(15) reads 



(p) 



if 



,(°) 



+ ie f (u { -y -u^ + ) =0 



(0) . . , ( J. ) 

— 7?i_l + ief[u_ + 



) = (36) 



-{-r/2 + i0)u { +y +ie f ( 



(o) 

u++ — u 



(o) \ _ 



)=o 



where 7 = Ea.b S ( e a)S(E b - 2/i)|(yl|7r| J B>| 2 and 6 = 
J2a.b S ( e a)\(A\tt\B)\ 2 2oj/{E 2 b -uj 2 ), which give the asymp- 
totic solution of the well-known optical Bloch equations 
in the rotating wave approximation. We recall that these 
equations are obtained by neglecting the non-secular terms 
of the Redfield part of (34) and the non-resonant cou- 
pling terms to the monochromatic field, which is valid in 
the limit of weak coupling to the heat bath and to the 
monochromatic field [26]. 



B Heat bath influence 

The functions jki.k'i' which appear in (11), are related to 
the time functions K k i tk ip given by (32), by 



lkl,k'l'{z) 



/ dte izt K kWl ,(t). (37) 
Jo 



(38) 



where 5(oj) = 5(lo) + (i/7r)u> -1 , since [f2, Hs] = 0. In this 
expression, E A and \A) denote the eigenvalues and eigen- 
states of the bath Hamiltonian H B , Pa = Tr[{2\A)(A\] is 
the initial population of state \A). For the initial state (3), 
P A = Z^cxpi-EA/T). 

For to = 0, V = k' and £/ = E k , expression (38) sim- 
plifies to 



Jkl,h 



k ,(i0+) = Y j Pa{-^ A MB){B\^ 1 \A) 

A.B j 



tt(<W + 6 kk >)6(E A - E B + uj kj ) + i- 



8lk' — $kk' 



Ea — E B + ujkj 
+ 2ir(A\ir k >i\B)(B\ir kk/ \A)6(E A - E B +Lo k , k )\, (39) 



which leads to 

Tki = lk k,u(^ + )e- E ' /T (40) 
= 2TTj2 p Ae- El/T \\(A\TTi k \B)\ 2 5(E A - E B +u lk ) 

A,B 

-5 kl Y,\{^ko\B)\ 2 5{E A -E B +u3 kj )^. 



For P A oc exp(-E A /T), since uj ik = E x - E k , T tk = T kx . 
For e = 0, ttu — en(k\S x \l) where S x = J2 n a n> and bence, 
for k ^ I, T k i oc S k i where S k i = (k\S x \l) 2 . Using (39), one 
finds, for Ei = E k , 

YjkLk'k'{iO + )pk' = -it P A (A\ir kk >\B)(B\ir k n\A) 

k< A.B.k' 



Pi - Pk 



Pk+Pi - Z-prPk' ) 5(u} A b +u kk >) + -- 

Pa J 7T was + uJkk' 



(41) 



where ujab — Ea — E B . For e = 0, the above summand 
vanishes if S kk < = 0, S k n = 0, or p k > = exp(uj kk > /T)p k = 
exp(uj kk > /T)pi, which ensures that (22) is the solution to 
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(20). Similarly, one finds, for Ey — E k > , 

J2Pk"fkk,k>l>(iQ + ) =7T PA(A\TTl'k\B)(B\TT kk ,\A) 

k A,B,k 

i Pk> - Pi' 



(2p k - p k , - pi,)5{u AB + cjfe/fe) + 



7T UlAB + Wfe'fe 



(42) 



For e = 0, this sum vanishes if p k = pi when Ski 7^ 0, 
which gives the left eigenvectors # g introduced in section 
5.5. 

The elements of the matrix Gq P \ defined in section 
5.5, are given by 



(p) 



= lki,k'v{iQ + + Uhi) 



= ire 2 p A I (A\n\B) \ 2 <2si n s kk >6(uj ABV 

A.B ^ 



) 



- ^2 \_^l' s kj s jk'S{ UJ ABkj) + S k k'Sl'jSjl5(—L0ABlj) 
3 



(43) 



where s k i = (k\S x \l) and uo AB kl = u AB + ^ki = E A - 
E B + E k — Ei. The equalities ui k i = ^k'v = P^ have been 
used to simplify this expression. The integer p does not 
appear explicitly in (43), but this expression is meaning- 
ful, for a given p, only for pairs (fc, I) and (k' , I') such that 2N 
ujki = w k >i> = puj. Since s k i = when k and I do not belong = 2 ^ A r Jr\, 

to the same set £ q , the element (43) is nonvanishing only r n =i 

when k and k' belong to the same set, and I and /' also. It 
can be seen, by writing in block diagonal form, that 
this property is also satisfied by the inverse matrix F^ p ' = 
(G P ^) _1 . Another useful property of the matrices Gq P ' is 
the following. We see that permuting k and I, and k' and 
V in (43), is equivalent to complex conjugation. Conse- 



is the dimension of the subsystem n Hilbert space [35,36], 
as follows. Since 

{$ r \W\$ r >) oc (V'i r) |^i r ' ) ) = 5 rr > where 
|*r> = ®£Lil^4 r) >, and similarly for H 1: A = (V\H\V) = 
E r l A r| 2 (^r|^|^r)- Consequently, A > (0\H\0), as this 
last value is the minimum possible one for the average 
energy, and |0) is unique. Thus, |0) is not 



D No decoherence-free subspace 

We show here that, even in the limiting case e/ = e = 0, 
there is no decoherence-free subspace in the TLS Hilbert 
space. This can be directly proved from the correspond- 
ing Hamiltonian expression. For tf = e = 0, the Hamil- 
tonian (1) simplifies to H — Htls + He + eirS x where 
S x = J2 n °ni an d hence, the coupling of the TLS to their 
environment, is described by a product term irS x where tt 
is an observable of the environment and S x is an observ- 
able of the TLS system. Thus, a decoherence-free subspace 
S would be a space spanned by eigenvectors of S x with 
the same eigenvalue and invariant under the TLS Hamil- 
tonian Htls- In such a vector space S, there would exist 
a basis consisting of eigenstates of Htls- In other words, 
Htls and S x would have common eigenvectors. Such a 
state is also eigenvector of [Htls, S x ] with eigenvalue 0. 

For any state \ip) = ^ r A r ®„ \Vn^)n, the state \ip') = 
[Htls,Sx]\iI>) is given by 



(r) 



v'/n-l 



"n+lj 



KvP] <\Vr) 
(44) 



Vn ^)n, with the conventions /i2jv+i-n 



quently, (G p) )ik.i'k' = (G p) )^ fc , ; , , and similarly for the 



inverses F^ p ^. At zero temperature, (G^)ki,k'i' vanishes 
if Ei > Ey, and hence (F^) k i,k'i' = for E t > E v . Thus, 
in particular, (FW) fei)C - a = 6 kc - Si a /(G^) C 



)c a,c a- 



where \rj r ) = 

(r) (r) 

—h n and r/_{ = VI^n+i = 0- We now assume that \ip) is 
an eigcnstatc of Htls- In this case, all the configuration 
states \rj r ) have the same energy. Since, for h n ^ 0, no 
coefficient in the decomposition (44) vanishes, a term in 
this sum does not contribute only if it is the opposite of 
another one. This is only possible if there exist n' and 
r' such that cr*,\r] r ,) — cr*\r} r ). For general values of the 
fields h n , \rj r ,) = (Jn'VnlVr) an d \V r ) can nave tne same 
energy only if n' = 2N + 1 - n, J^jv+i-n = Vn \ and 



(r) 



(r) 

7 l2N-n+2 



.JT) 
'ln-1 



C No thermal equilibrium Schmidt 
decomposable state 

We consider a composite system consisting of Af subsys- 
tems, with no long-range interaction between these sub- 
systems. The Hamiltonian H of the complete system can 
thus be decomposed as % = "Hi + W where no observable 
of subsystem 1 appears in %' and there exists a subsys- 
tem n^l which is not affected by H\. The Hamiltonian 
H is assumed to have a nondegenerate ground state |0). 
The thermal equilibrium state of the sytem is thus pure 
in the zero temperature limit. We show that |0) is not a 

Schmidt decomposable state — Ylt=i^r ®n=i IV^) 

where (ipn^ipn = cW, s g {2, . . . , min„(d„)}, and d n 



<U if " i {N,N + 1} 

(r) (r) 

replaced by rj N+2 = —rj N _ 1 if n = N or N + 1. Con- 
sequently, for given n and r, if n! and r' exist, they are 
unique. Consider now a given r. The conditions on f] r 
found above, clearly show that there cannot exist n_i, 
r_i, no, ro, n\ and n such that (Jn T \ r lr T ) = ^n+rl^r)- 
Thus, there always remains at least one non- vanishing 
((Vr\Oi H TLS,S x }\iP). Therefore, [H TLS , S*]|V>> ¥= 0, and 
hence Htls and S x have no common eigenvector. There 
is thus no decoherence-free subspace in the TLS Hilbert 
space. This can happen in particular cases. For exam- 
ple, for N = 2 and h 2 = 0, (| + - + +) - | - - + -) - 

I + H h> -h | — I ))/2 is eigenstate of both H T ls 

and S x , and hence of H = Htls + H& + eirS x . If the TLS 
are initially prepared in this state, they remain in it for 
ever. 
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E Two-heat-reservoir environment 

In this Appendix, we show that the TLS asymptotic state 
Pco can be a pure Schrodinger cat state if the monochro- 
matic field is replaced by a second heat bath B' . We con- 
sider the Hamiltonian 



27V 



2N 



H = H TLS + H B + H B ,+eJ2a^ + e'J2<^n (45) 



where it is an observable of bath B and ir n are observablcs 
of bath B' and the initial state 



]r k i\k)(l\® Z 



-l e -H B /T 



-l e -H B ,/T' 



(46) 



k.i 



where Z' = Tr cxp(— H^'/T) and T' is the temperature of 
bath B'. We are concernerd with the regime e <C e' <C e, 
which is the analog of that considered in section 5.5. In this 
regime, the deviation of the coupling to bath B from the 
ideal form S x tt, can be neglected, and e has been set to zero 
in the above expression of H. Here, p,^ is a steady state. 
From section 5.2, we know that it is of the form (22). As in 
section 5.4, it can be shown that the potential probabilities 
p q where q ^ ±, vanish in the low temperature limit, faster 
than exp(— 3J/T). The difference is simply that here the 
downward transitions are induced by bath B' . 

The ratio p + jp_ can be determined by a perturbative 
calculation similar to that done in section 5.5. The state 
Poo obeys J2k>,i> Iki.k'V {i^ + ) u k% = 0, where k and I are 
such that Ek = Ei, k' and V are such that E^ = Ey , 
and r yki,k'i , ('iO + ) is given by an expression of the form (38) 
but with sums over eigenstates of both baths B and B' . 
This equation set can be written in matrix form as je 2 G + 

(e') 2 G']u = where G is the matrix Gq ^ introduced in 
section 5.5. The matrix term proportional to ee', vanishes 
under the assumptions Tr(.!?7r) = Tr(J77r^) = 0. Using the 
left and right eigenvectors of G, <f> q and ^Sf q , it can be 
shown that ^ , r qq >p q i = 0. Here the rates r qq i are given 
by, for q^q', 



2tt \ ^ x ^ 

q ke£ q A,B 



-E A /T'-E,/T 



e-^'HiEA-EB+uik) 



J2(l\<\k)(A\n*\B) , (47) 



where A and B run over the eigenstates of bath B', and by 



iq'q' 



E 



qjtq> T 11 



All transitions, induced by bath £>', 



from a state |Z) given by £ q > to a state \k) given by £ q , con- 
tribute to (47). For T' = T, z q >r qq i = z q r q / q , which ensures 
that the equilibrium density matrix p eq oc cxp(— Htls/T) 
is a steady state of the TLS system. We consider the low 
temperature regime T" -C T -C h n < J. In the sum 
(47), the terms such that Ei < E k , vanish in the zero 
T" limit, whereas the terms such that Ei > Ek, reach fi- 
nite values. Thus, in the regime considered, the upward 
transitions induced by bath B' are negligible and only the 



downward transitions contribute to (47). Consequently, 
r+_ - cxp[-(£ c - E )/T], r ++ ~ -r_+ < r+_ (for 
N > 3), and hence, since p q <€. exp(— 3J/T) for g^±,p_ 
vanishes, and the TLS steady state is |Scs + ). 
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